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Executive  Summary 


Building  on  our  previous  work,  the  main  objectives  of  this  research  were  to  further  develop  and  expand  the 
necessary  fundamental  algorit  hms  mid  procedures  to  enable  high-fidelity  multidisciplinary  design  of  complete 
aircraft  con  figurations.  Based  on  a  number  of  achievements  in  our  previous  work,  we  believe  that  wo  have  a 
mathematical  formulation  for  the  tools  that  are  required  to  create  such  an  integrated  design  environment. 
We  were  able  to  increase  the  scope  of  our  research  by  leveraging  from  other  sources  of  support.  In  particular 
Georg  May’s  thesis  work  benefited  from  the  support  of  a  Stanford  Graduate  Fellowship. 

In  order  to  realize  a  truly  automated  mult  ill  iscipltthtry  design  environment,  several  tasks  needed  to  be 
accomplished  to  enhance  the  tools  that  we  already  had  and  to  integrate  them  in  such  a  way  that  coordination 
between  the  evolution  of  the  disciplines  is  maintained. 

This  report  describ.es  the  progress  that  has  been  made  towards  tin*  accomplishment,  of  these  goals.  The 
fundamental  efforts  conducted  have  Ionised  on  four  different  areas,  Firstly,  we  have  continued  to  develop  our 
analysis  and  design  unstructured  mesh  capabilities  for  viscous  flows;  a  later  sect  ion  provides  a  summary  of  I  his 
work.  Secondly,  we  have  further  developed  our  viscous  planfonu  optimization  capability  which  now  includes 
a  model  of  the  wing  weight  to  guide  the  selection  of  the  planfonu  variables  (area,  aspect  ratio,  sweep*  taper 
ratio,  etc,)  Third,  we  have  carried  out  initial  efforts  to  c  reate  truly  high-fidelity  multidisciplinary  designs 
by  creating  a  two-level,  multi-fidelity  design  strategy  that  has  been  demonstrated  using  a  small  supersonic 
jet  as  an  example.  Finally,  we  have  explored  the  possibility  of  developing  au tomcat ic  feedback  control  of 
aerodynamic  flows  based  on  a  combination  of  adjoint  based  sensitivity  analysis,  reduced  order  modelling,  and 
linear  quadratic  control  theory.  A  procedure  for  controlling  rite  surface  pressure  disturbances  via  blowing  and 
suction  lias  been  demonstrated  and  this  has  been  extended  to  show  the  possibility  of  active  flutter  control 

1  he  theses  and  main  publications  which  have  resulted  from  the  research  performed  in  association  with 
the  grant  are  listed  below. 
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Unstructured  Methods  for  Analysis 
and  Design  of  Viscous  Flows 


In  this  section  we  report  our  progress  towards  a  new  platform  for  computational  aerodynamics  analysis 
and  design  on  arbitrary  meshes.  The  approach  was  designed  for  maximum  flexibility  to  serve  as  the  basis 
for  a  future  industrial  strength  How  solver  on  general  grids,  and  as  a  framework  for  advanced  research  in 
the  area  or  CFD  and  aerodynamic  design.  Such  a  flexible  platform  is  crucial  for  transfering  new  research 
to  industrial  applications  We  briefly  describe  the  capabilities  and  methods  and  show  results  from  initial 
validation  on  in  viscid  and  viscous  test  cases.  The  research  addresses  open  issues  in  CFD  on  unstructured 
grids  such  as  viscous  discretization  and  multigrid  methods  in  an  unstructured  context. 


Support  For  Arbitrary  Meshes  and  Different  Discretization  Techniques 


Figure  I:  Examples  of  different  mesh  types  mid  discretization  techniques.  From  left  la  right:  Structured, 
unstructured  in  cell-centered  discretization  (variables  are  stored  at  the  cell-centers),  unstructured  in  cell- 
vertex  discretization  using  the  median  dual  { variables  are  stored  at  the  nodes),  Cartesian, 

Different  types  of  meshes  (structured,  unstructured,  Cartesian)  as  well  as  conceptually  different  methods 
of  discretization  are  used  in  computational  fluid  dynamics  today.  Examples  of  discretization  techniques  are 
cell -centered  discretization  or  cell- vertex  discretization  wiili  one  of  the  possible  ways  i  n  define  dual  meshes 
from  a  given  primary  tnesli  topology.  Figure  i  illustrates  different  mesh  types  and  different  discretization 
techniques. 


Figure  2:  Two  dimensional  illustration  of  the  face-based  data  structure. 
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do  i=ncsill  Pncell2 

set  reaidualti)  to  zero 
end  do 

do  n*nf acel pnfaee2 
K\  -  ncUltn) 

N2  -  ncf  <2 , n) 

flux(N)  *  f (solutionCNl) Psolution(N233 

residual CN2)  =  residual(N2)  +flux(N) 
residual (Nl)  =  residual (Nl)  -flux(N) 
end  do 

do  i-ncelU  ,nceU2 

solution Ci)  *  solution(i)  -residual f i) 
end  do 


Figure  1  :  Stencil  for  2nd  order  flux  compulation  oil  a  structured  mesh. 


The  newly  developed  method  can  operate  on  any  of  the  above  mesh  types  and  discretization  technic  pies. 
In  fact,  the  flow  solver  module  is  at  present  completely  indifferent  to  the  underlying  mesh  topology  so  that 
it  can  rum  in  principle,  on  arbitrary  meshes.  A  face -based  data  structure  is  used,  which  recognizes  pointers 
from  cell  interfaces  to  adjacent  nodes  as  the  only  connectivity  information.  This  is  illustrated  in  Figure  2. 
Regardless  of  the  mesh  topology,  there  will  always  be  a  finite  number  of  interfaces  separating  exactly  two 
control  volumes,  including  halo  cells  at  boundaries.  By  looping  over  these  interfaces  all  necessary  flux 
computations  can  be  carried  out.  The  pseudo-code  in  Figure  illustrates  the  basic  algorithm  used  in  the 
How  solver  module.  While  this  algorithm  is  simple,  the  connectivity  information  used  is  not  normally  the 
one  stored  in  primary  meshes.  Computational  complexity  is  shifted  to  the  preprocessing  stage,  where  the 
necessary  data  structure  is  generated  along  with  metric  data.  Currently,  preprocessing  is  implemented  to 
provide  metrics  for  the  median  dual  mesh  and  the  primary  mesh  for  the  four  most  common  element  types 
tetrahedrm  prisms,  pyramids,  and  hexahedra. 

Preprocessing  for  hexahedra  allows  the  treatment  of  structured  meshes  in  an  unstructured  context.  For 
structured  meshes  the  data  structure  can  also  be  augmented  to  recover  the  treatment  that  is  normally  used 
For  this  type4  of  mesh.  Figure  1  illustrates  the  usual  stencil  used  for  a  second  order  flux  compulation  on  a 
structured  mesh.  The  blue  cell  centers  in  Figure  1  are  pointees  of  the  face  between  them  in  an  unstructured 
context.  By  creating  an  additional  level  of  pointers  to  the  next-but-ono  neighbors  (green  ceil  centers),  it  is 
possible  to  reproduce  the  structured  algorithms  with  the  unstructured  code.  Note  that  there  me  thus  two 
ways  in  which  structured  meshes  can  be  treated.  One  way  is  to  look  at  them  as  unstructured  hexahedrnl 
meshes  and  use  a  general-mesh  treatment.  Involving  gradient  reconstruction  (to  be  described  below).  Another 
way  is  to  use  the  described  extra  data  structure  to  reproduce  the  exact  structured  algorithms,  which  is  faster 
compared  to  lit  si  method.  It  also  serves  to  asses  the  minimum  overhead  that  is  incurred  compared  to  a 
structured  code  due  to  the  indirect  addressing.  Numerical  experiments  have  shown  this  overhead  to  be 
around  25%. 

Cell-centered  discretization  using  the  primary  mesh  or  cell-vertex  discretization  using  the  median  dual 
mesh  can  be  selected  at  run-time,  ihe  only  difference  being  the  preparation  of  the  metric  data.  Usually, 
subtle  differences  between  cell -centered  and  cell-vertex  discretization  exist  even  in  a  mesh  transparent  flow 
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Figure  5:  Pressure  distribution  in  the  wing  section  y/b  =  0,2  for  the  Oueni  Mti  wing  at  Mach  number 
Af  —  (1,81  and  a  —  3,06"  angle  of  attack.  Results  art*  shown  fur  a  structured  mesh  with  1.1  million  nodes 
ami  for  a  tetrahedral  mesh  with  316.000  nodes, 

solver  module,  for  example  at  t  he  boundaries.  In  this  implementation  a  unified  treatment  of  the  boundaries 
has  been  chosen,  using  halo-cells  at  boundaries  for  any  kind  of  mesh.  i.e.  for  cell-vertex  schemes  boundary 
conditions  are  not  enforced  at  the  nodes.  During  preprocessing  boundary  nodes  become  dual  cell  centers, 
which  are  no  longer  located  on  the  boundary,  but  shifted  to  the  dual  cell  centroids.  The  construction  of  lmlo 
cells  proceeds  in  the  same  manner  as  for  cell-centered  discretization. 


In  viscid  Validation 


Figure  5  provides  initial  validation  of  the  uni  fit?d-  mesh  algorithm.  Both  a  structured  mesh  with  approximately 
11  million  nodes  and  an  unstructured  tetrahedral  mesh  with  316,000  nodes  have  been  used  to  calculate  the 
How  for  the  Onera  M6  test  ease  at  Mach  number  A!  ~  0,81  and  q  =  3.00"  angle  of  attack.  The  cp  distribution 
at  the  wing  section  y/b  —  0.2  is  compared  to  results  obtained  with  the  solver  AIRPLANE  on  the  same 
unstructured  tetrahedral  mesh.  All  computations  have  been  carried  out  using  cell-vertex  discretization  on  a 
three-level  nniltigrid  sequence. 


Viscous  Discretization 


Since  gradients  are  needed  for  state  reconstruction  they  are  initially  computed  for  the  nodes  of  the  primary 
mesh  or  the  centroids  of  the  primary  elements,  depending  on  the  discretization  type  chosen.  For  viscous 
dieretization  the  gradients  must  be  evaluated  at  the  interfaces  between  adjacent  cells.  This  can  be  dune  by 
u  suitable  average,  which  typically  depends  on  the  mesh  structure  and  choice  of  control  volume.  Since  ihe 
goal  for  the  computational  platform  is  complete  mesh  transparency,  tins  dependence  on  the  mesh  topology 
is  unsatisfactory. 

In  t lie  current  implementation  gradients  are  constructed  using  either  the  unweighted  least -squares 
method  or  a  gauss-type  reconstruction,  i.e. 


U) 


where  \\  is  the  cell  volume,  and  tin*  overliar  denotes  a  suitable  face  average. 


For  conventional  viscous  discretization  tin  initial  value  for  the  gradient  at  the  cell  interface  may  he 
obtained  by  averaging  the  values  of  adjacent  cells,  o,  and  6,.  to  give 
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Subsequently  a  correction  in  the  direction  is  applied  to  avoid  odd-even  decoupling,  where 
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This  leads  to  the  following  gradient  on  the  interface 
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This  modilicatiou  results  in  a  stable  diseret i/al ion  of  the  viscous  terms  and  does  not  ali'eet  the  accuracy 
its  long  as  the  interface  integration  point  is  located  halfway  between  the  two  cell-centers.  This,  however,  is 
only  tin1  ease  for  the  median  dual  control  volume,  rendering  this  type  of  viscous  discretization  highly  mesh 
dependent. 


Using  concepts  of  kinetic  gas  theory,  viscous  discretization  ran  be  accomplished  in  a  different  fashion 
Many  researchers  lmvc  used  the  Boltzmann  equation  or  its  BGK  sinipHfieatioti,  he.  the  BGK  equation,  as 
a  basis  Tor  the  construction  of  numerical  fluxes*  Firstly  we  note  that  the  Navier-Stokes  equations  can  he 
obtained  by  Chaptnan-Enskog  expansion  of  both  the  full  Boltzmann  equation  as  well  as  the  BGK  equation* 
the  difference  being  only  in  the  values  of  the  transport  coefficients. 

For  more  details  of  the  actual  discretization  used,  the  reader  is  referred  to  the  appropriate  publication 
at  the  end  of  this  report .  Tire  details  have  been  omitted  here  for  conciseness. 

For  initial  validation  we  provide  the  solution  of  a  resolved  shock  structure  obtained  with  the  new  scheme, 
see  Figure  tb  It  should  tie  pointed  out  at  this  point  that,  since  the  Chapmami-Btiskog  expansion  of  the  BGK 


Shoe*  Sfructuri  CompulBlion 


Velocity  distribution  across  a  normal  shock. 


Heat  flux  across  a  normal  shock. 

Figure  (>:  Resolved  shock  structure  for  upstream  Mach  number  A/  =  1  At  and  a  PrandtJ  number  of  Pr  —  0.72 
equation  gives  the  transport  coefficients  as  a  sole  function  of  the  collision  time  r,  only  one  parameter  can 


he  set  automatically  to  the  right,  value.  This  will  usually  he  the  dynamic  viscosity.  For  accurate  heal  flux 
computations,  the  reconstruction  involving  the  temperature  gradient  lias  to  be  scaled  appropriately,  so  as  to 
correspond  to  the  correct  Praiidtl  number. 

For  further  validation  consider  a  zero- pressure-grad  lent  boundary  layer.  Figures  7  shows  that  the  dilu¬ 
tion  computed  is  in  close  agreement  with  the  Bhssius  solution, 

A  detailed  description  of  this  research  is  given  in  the  thesis  of  Georg  May,  entitled  "A  Kinetic  Scheme 
for  the  Navicr- Stokes  Equations  and  High-Order  Methods  for  Hyperbolic  Conservation  Laws".  During  the 
period  of  Ins  thesis  research  Georg  May  received  support  from  a  Stanford  Graduate  Fellowship. 
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Figure  7:  Laminar  Boundary  Layer  at  lie  —  GOOD  and  M  =  1)  2  using  the  1ST  scheme  and  direction-corrected 
central  discretization  of  the  viscous  terms. 
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Planform  Optimization  Studies  Using 
Viscous  Flows  and  a  Weight 
Estimation  Method 


In  this  section  we  focus  on  our  work  on  aero-  struct  oral  optimization  of  wings  for  long  range  transport 
aircraft,  using  adjoint-based  optimization  techniques.  We  arc  interested  in  exploring  and  comparing  the 
attainable  irado-offs  such  as  L/l)  vs.  Mach  number,  which  may  be  appreciable  higher  than  the  historical 
trends  typically  used  in  conceptual  design.  We  also  seeks  to  identify  a  discernable  trend  in  the  variation 
of  planform  variables  such  as  sweep,  thickness-to- chord  ratio,  aspect  ratio,  and  chords  for  optimum  wings. 
Results  form  wing- fusel  age  and  eo  m  pie  1e-aircr.MWo»  figuration  optimizations  indicate  that  stretching  the 
span  together  with  decreasing  the  sweep  and  thickening  the  wing  sections,  the  lifMu-drag  ratio  can  be 
increased  without  any  penalty  on  the  structure  weight. 

J  lie  methodology  used  for  this  wHght-estimate-based  optimization  studies  has  been  presented  earlier 
and  will  not  be  repeated  here.  The  reader  is  referred  to  the  publications  by  Leoviriyakit  and  Jameson  listed 
at  the  end  of  this  report. 

Redesign  of  the  Boeing  717  wing 

We  present  a  result  to  show  that  the  optimization  approach  we  are  pursuing  can  successfully  trade  plan  form 
parameters  and  wing  weight.  We  also  demonstrate  how  to  apply  strategy  game  theory  to  gradient  based 
optimization. 

Here,  the  ease  chosen  is  the  Boeing  747  wing  fuselage  combi  nation  at  Mach  U.87  and  a  lift  coefficient 
Cl  —  0.42.  The  computational  mesh  is  shown  in  Figure  8. 

[n  f  his  test  case,  the  Mach  number  is  the  current  normal  cruising  Mach  number  of  0.85.  We  allowed 
section  changes  together  with  variations  of  sweep  angle,  span  length,  chords,  and  section  thickness.  Figure 
9(a)  shows  the  baseline  wing.  Figure  9(b)  shows  the  redesigned  wing.  The  parameter  ^  was  chosen  such 
that  the  cost  function  corresponds  to  maximizing  the  range  of  the  aircraft.  Here  in  30  design  iterations  the 
drag  was  reduced  from  137  counts  to  117  counts  and  the  structural  weight  was  reduced  from  198  counts 
{80.480  Iks)  to  HM  counts  (75,000  lbs).  The  large  reduction  in  drag  is  the  result  of  the  increase  in  span  from 
212. 1  ft  to  231.7  ft,  which  reduces  the  induced  drag.  The  redesigned  geometry  also  has  a  lower  sweep  angle 
and  a  thicker  wing  section  in  the  inboard  part,  of  the  wing,  which  both  reduce  the  structural  weight.  Moreover 
the  section  modification  prevents  i  lie  formation  of  shock.  The  baseline  and  optimized  planforms  are  shown 
in  Figure  9(c).  Overall,  the  re-design  with  variation  planform  gives  improvements  in  both  aerodynamic 
performance  and  structural  weight,  compared  to  the  previous  optimization  with  a  fixed  planform. 


Redesign  of  the  BAe  MDO  DATUM  wing 

To  further  validate  this  planform-and-sectibn  trend,  we  selected  the  BAe  MDO  DATUM  wing,  a  wing  that 
has  been  designed  with  relatively  modern  tools.  At  its  cruising  Mach  .85,  this  wing  has  low  sweep  angle  and 
high  thickness- to-ehord  ratio  sections. 

This  test  case  presents  a  technical  challenge  to  the  optimization  because  the  BAe  and  B747  are  designed 
to  operate  at  the  same  flight  condition  and  their  planforms  are  sized  in  the  same  range.  However  the  original 
sweep  of  BAe  is  already  smaller  than  the  optimum  sweep  of  B7  17  and  its  wing  span  is  already  longer  than 
the  optimum  span  of  B747. 

Figures  10(a),  10(b),  and  10(c)  show  the  original  wing,  optimized  wing,  and  their  planforms  respectively 
Despite  of  low-sweep,  long-span,  and  thick- wing-sec  turns  of  the  original  wing,  the  optimal  wing  lias  less  sweep. 
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Figure  8:  Computational  grid  of  tin1  13717  wing  fuselage 


longer  span,  and  thicker  wing  sections.  Bur  tin*  changes  in  the  planform  are  not  large.  With  these  changes, 
the  optimum  wing  shows  improvement  in  both  drag  and  weight.  The  drag  is  reduced  from  H>1  counts  to  I  15 
counts,  and  the  weight  is  reduced  from  480  counts  (87.500  lbs)  to  170  counts  (86,980  lbs).  This  optimized 
BAe  wing  strongly  agrees  with  the  trend  suggested  from  the  B717  and  MD11  cases. 

This  research  is  documented  in  t  he  thesis  of  Kasidit  Leoviriyakit.  entitled  “Wing  Planform  Optimization 
via  an  Adjoint  Method", 
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Figure  S);  Redesign  of  Boeing  747  wing*  using  section  and  plan  form  modifications.  The  optimum  wing  has 
longer  spun*  loss  sweep,  and  thicker  wing  sections.  We  also  over-plot  the  optimum  plan  form  form  our  in  viscid 
culru  hit  ion  to  indicate  good  agreements  between  tin1  in  viscid  ami  viscous  optimizations.  Top  left:  baseline. 
Top  right;  Optimized,  Plan  form  {Baseline  :  Green.  Optimum  :  Bine). 
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Figure  10:  Redesign  of  RAe  MDO  DATUM  wing,  using  section  and  plnnlbrm  modi  Heat  ions.  The 
wing  has  longer  span,  less  sweep,  and  thicker  wing  sections. 
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Two-Level,  Multi-Fidelity  Approach 
to  Multi-Disciplinary  Design  of 
Supersonic  Aircraft 


The  ronatptual/pndhninajt/  design  of  supersonic  aircraft  configurations  requires  a  multi- disci  pi  man 
approach  that  provides  the  designer  with  information  regarding  the  key  trade-oils  between  the  disciplines 
participating  iti  the  design.  At  the  same  time,  at  these  stages  of  the  design,  the  available  tools  must  provide 
a  level  of  flexibility  i  hat  permits  i  he  exploration  of  large  areas  of  Hie  design  space  with  significant  changes  lo 
a  baseline  configuration.  In  order  to  achieve  credible  results  one  would  like  to  use  high-fidelity  modeling  tools 
for  all  of  the  components  (and  interactions)  of  the  design.  This  can,  however,  be  prohibitively  expensive  and 
in  addition,  it  may  significantly  decrease  the  ability  lo  make  drastic  modifications  to  the  aircraft  configuration 
in  question .  As  our  work  has  progressed  in  this  area,  we  have  come  to  realize  that  a  truly  hybrid,  multi- 
fidelity  approach  that  is  properly  managed  is  one  of  I  he  answers  to  the  supersonic  design  problem.  We  are 
currently  pursuing  a  two-level  approach  to  the  design  of  a  supersonic  business  jet  configuration  where  we 
combine  a  conceptual,  SIMPLEX-lmsed,  low-fidelity  optimization  tool  with  a  hierarchy  of  flow  solvers  of 
increasing  fidelity  (including  simplified  aerodynamic  models,  a  linearized  panel  method  and  both  structured 
and  unstructured  Euler  solvers)  and  advanced  adjoint-based  Sequential  Quadratic  Programming  (SQP) 
optimization  approaches.  Although  this  kind  of  aircraft  has  been  studied  in  the  past  within  the  context 
of  low  supersonic  boom,  in  this  work  we  focus  on  the  aerodynamic  performance  aspects  alone;  no  attempt 
is  made  to  reduce  the  acoustic  signature  since  this  has  tad  been  a  driver  in  military  applications  of  this 
technology.  The  results  show  that  this  particular  combination  of  modeling  and  design  techniques  is  quite 
effective  to  produce  designs  with  optimum  performance  that  meet  or  exceed  all  of  the  design  constraints 
of  l he  problem.  In  addition,  wc  show  that  high-fidelity  aerodynamic-  shape  optimization  techniques  for 
complex  configurations  (such  as  the  adjoint  method)  can  be  effectively  used  within  the  context  of  a  truly 
1 1  mlti-discipli  i  wry  design  eu vi  run  r non t . 

In  i  his  work,  wc  combine  idea*  of  mull  i-fidelhy  analysis  and  design  and  a  two- level  oplhniz&tkm  proce¬ 
dure  into  a  hybrid  concept  that  includes: 

1  The  Program  for  Aircraft  Synthesis  Studies  (PASS):  n  multi-disciplinary  design  tool  that  incorporates 
carefully  tuned  fast  models  for  the  various  disciplines  in  the  design  and  is  able  to  deal  with  all  the 
major  objective  functions  and  constraints  in  typical  aircraft  synthesis  problems. 

2.  A  hierarchical,  multi-fidelity  response  surface  generation  technique  that  uses  results  from  classical 
supersonic  aerodynamics,  a  linearized  supersonic  panel  code  ( A 502/ Panair),  and  unstructured  adaptive 
Euler  solver  ( AirplanePhis)  to  create  models  of  the  aerodynamic'  performance. 

il.  Automated  tools  based  on  a  coratmm  geometry  database  to  drive  the  analysis  tools  that  are  used  in 
the  generation  of  the  response  surfaces  in  this  problem  (BOOM-UA),  This  CAD-tu-solutkm  procedure 
is  based  on  the  CAPRI  CAD- interface  of  Haimes,  the  A50'2/Panair  and  Airplane  PI  us  flow  solvers,  mid 
the  Cental ir  mesh  generation  system, 

1.  Adjoint  aerodynamic  shape  optimization  tools  for  both  single-block  wing-body  configurations  (SYX87- 
SB)  and  mult  iblock  complete  configurations  (SYN107-MB)  that  use  inexpensive  gradient  calculations 
with  larger  numbers  of  design  variables  to  modify  the  twist  and  camber  distributions  of  the  wing 
(without  changes  to  the  wing  plnnform)  and  to  achieve  the  highest  aerodynamic  performance. 

More  details  of  each  of  the  components  of  this  work  can  be  found  in  the  appropriate  publications  listed  at 
t he  end  of  this  report, 

A  typical  unstructured  surface  mesh  for  the  types  of  aircraft  configurations  in  question  t  an  be  seen  in 
Figure  I  1  below.  I'he  resulting  flow  solution  is  shown  in  Figure  12. 


Figure  11;  Unstructured  tetrahedral  surface  mesh  around  full  baseline  configuration 


Multi-fidelity*  Multi-Discipline  Design  Approach 

This  section  explains  the  procedure  we  have  used  to  integrate  the  tools  described  above  into  a  single  analysis 
and  optimization  capability.  The  concept  is  straight  forward:  if  the  multi-fidelity  analysis  capability  can 
he  used  to  create  response  surfaces  for  the  drag  coefficient,  C'q T  the  corresponding  low-fidelity  modules 
in  PASS  can  be  replaced  by  these  response  surface  fits,  [’his  makes  for  a  remarkably  simple  integration 
problem  and  also  provides  us  with  the  ability  to  predict  the  changes  in  aerodynamic  performance  resulting 
from  wing  section  changes.  The  baseline  version  of  PASS  is  unaware  of  the  actual  wing  sec  tions  used  and 
assumes  that,  whatever  the  sections  arc.  they  have  been  adjusted  in  such  a  way  that  the  camber  and  twist 
distributions  are  optimal  (in  the  sense  that  they  get  dose  to  elliptic  load  distributions  in  both  the  span  wise 
and  slreamwise  directions).  PASS  can  then  be  used  to  generate  optimized  results  and  the  outcome  of  the 
optimization  can  be  analyzed  using  the  high-fidelity  tools  to  ensure  that  the  response  surface  lits  provide 
accurate  representations  of  l lie  true  high-fidelity  responses.  The  level  of  accuracy  iu  the  response  surface 
representation  depends  greatly  on  the  number  of  high-fidelity  calculations  that  are  used.  Since  we  are  trying 
to  minimize  this  number,  we  will  undoubtedly  incur  some  errors  in  the  fits.  The  validity  of  these  fits  is 
typically  assessed  by  direct  analysis  of  the  resulting  optimized  designs  using  two  different  Ruler  solvers, 

Oju  multi-fidelity  approach  to  the  construction  of  t  he  response  surface  fits  relies  on  a  hierarchy  of  three 
different  aerodynamic  analysis  modules 

L  PASS  internal  analysis  based  on  classical  aerodynamics. 

2.  A  502/Panair  supersonic  linearized  panel  code. 

3.  fouler  solutions  of  the  highest  fidelity  using  unstructured  mesh  (with  a  total  of  around  U2  million  nodes 
lor  the  complete  configuration)  -  we  refer  to  these  compulations  by  the  label  "Fine  Euler  (FE)." 

In  order  lo  obtain  response  surface  fils  of  the  highest  fidelity  one  could  carry  out  a  large  number  id  FE 
solutions  and  fit  the  resulting  data.  Unfortunately,  for  large  dimensional  design  spaces  (we  will  he  using  a 
total  of  23  design  variables  later  on),  accurate  fits  require  a  large  number  of  function  evaluations.  This  is 
particularly  true  in  our  case  since  the  ranges  of  variation  of  each  of  the  design  variables  will  be  rather  large. 
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Figure  12:  Upper  surface  pressure  distribution  for  full  baseline  configuration  using  AirpliinePlus  Ruler  cal¬ 
culation. 

The  main  objective  in  this  section  is  to  generate  response  surface?  lit  s  of  the  same  quahty /accuracy  i  hat 
would  be  obtained  by  evaluating  the  FE  solutions  only,  but  at  a  ranch  reduced  cost.  We  accomplish  this  by 
relying  on  a  fundamental  hypothesis  that  will  be  tested  later  on:  the  higher  fidelity  tools  art  only  nettled  tn 
small  region#  of  the  design  space  where  the  tower  fidelity  model#  have  exhausted  their  range  of  applicability. 
This  is  bound  to  lie  true  as  it  is  the  premise  upon  which  aerodynamic  design  has  been  predicated  for  the  last 
50  years:  aerodynamicists  and  engineers  use  the  fastest  tools  for  a  specific  purpose  (when  they  are  known 
to  work  well)  and  switch  to  more  time-consuming^  expensive  tools  only  when  they  are  needed.  For  example, 
in  supersonic  design*  classical  equivalent  area  concepts  and  linearized  panel  codes  can  provide  very  accurate 
results  as  long  as  non-linear  efl'eels  (such  as  transonic  flows  in  the  direction  normal  to  the  leading  edge  of 
the  wing)  are  not.  present  and  viscosity  does  not  play  a  dominant  role  in  i  lie  solution  of  the  flow* 

Willi  this  in  mind,  wo  have  used  the  following  five-step  procedure  to  create  the  response  surfaces  used 
in  this  work.  All  databases  of  candidate  designs  are  obtained  by  populating  the  design  space  using  a  Latin 
Hypercube  Sampling  (LHS)  technique, 

L  Run  a  large  database  of  candidate  designs  (>  8,000)  using  the  aerodynamics  module  in  PASS.  Each 
evaluation  takes  roughly  1  second  to  compute  on  a  modern  workstation  (Pentium  4.  3,2  GHz),  This 
evaluation  also  flies  each  Aircraft  through  the  mission  and  returns  a  measure  of  the  iufeasibilitv  of 
the  design  (an  L  2  norm  of  the  constraint  violations,)  Those  designs  that  are  found  to  significant b 
violate  the  requirements/constraints  of  the  mission  are  removed  from  the  database  and  are  no  longer 
considered  in  the  response  surface  creation. 

2,  Run  the  remaining  database  of  candidate  designs  2.  500)  using  the  A502/Panair  solver.  Each 
evaluation  requires  about  10  seconds  of  CPU  time  on  the  same  modern  workstation, 

3.  Select  the  design  points  whose  relative  error  for  ('n  (based  on  the  baseline  design)  is  larger  than  a 
specified  t  hreshold,  t  r.\$s-A502y  and  analyze  only  those  designs  using  the  Fine  Euler  (FE)  approach. 
In  our  work,  we  have  set  this  threshold  to  about  15  %  resulting  in  a  number  of  high-fidelity  function 
evaluations  in  the  neighborhood  of  200.  Each  FE  evaluation,  from  beginning  to  end,  including  geometry 
and  mesh  generating  the  bottlenecks  in  the  process*  since  they'  are  run  serially}  requires  about  10 
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minutes  oi  wall  clock  time,  The  flow  solution  portions  (using  AirphmePlus)  are  run  in  parallel  using 
Hi  Athlon  AMD210D+  processors  of  a  Linux  Beowulf  duster. 

1  A  baseline  quadratic  response  surface  fii  (using  least  squares  regression)  is  created  for  the  ( 'u  obtained 
with  A502/Pftiiair,  The  error  between  the  values  of  the  FE  evaluations  and  the  predictions  of  these 
quadratic  fits  is  approximated  with  a  Knging  method,  and  the  resulting  approximation  is  added  to  the 
baseline  quadratic  ills. 

In  sunn  the  response  surfaces  provided  to  PASS  are  the  addition  of  t  he  quadratic  fits  based  on  t  he* 
A502/Panair  results  and  the  Knging  fits  of  the  error  between  the  FE  solutions  arid  those  quadratic  fits. 

Figure  13  shows  the  result  of  the  over  2,500  candidate  designs  (green  dots)  evaluated  using  A502/Pauair 
that  are  retained  after  the  initial  filtering  of  over  8,000  PASS  results.  The  red  dots  in  the  Figure  indicate 
those  candidate  designs  for  which  the  predicted  values  of  C&  are  ofT  by  more  than  tpASS  A?> m  >  15% 
between  PASS  and  A 502,  Note  that  a  number  of  these  red  dots  have  unreasonably  large  values  of  CD  since 
the  geometries  and  design  conditions  are  such  that  the  limits  of  applicability  of  A502  are  exceeded,  These 
points  for  which  the  disagreement  between  PASS  and  A5CJ2  is  large  are  taken  for  further  evaluation  using 
FE,  Figure  14  shows  in  blue  the  results  of  the  FE  analyses  for  a  subset  of  about  200  of  the  red  A502  results. 
Hie  final  result  is  a  set  of  VK  evaluations  that  are  meant  to  be  clustered  around  the  areas  where  the  lower 
fidelitv  models  cannot  accurately  predict  the  flow  physics. 
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Figure  }‘A:  Database  of  PASS  (green)  and  A  502/Panair  (red)  results. 

This  inult  i-fideliiy  procedure  has,  to  some  extent .  the  flavor-  of  Richardson's  extrapolation  in  that  it 
recursively  uses  results  from  different  fidelities  to  arrive  at  a  final  answer /fit.  It  also  has  an  adaptive  nature 
to  it,  as  results  from  t  he  higher  fidelity  models  are  only  evaluated  in  areas  of  the  design  space  where  the  lower 
fidelity  models  are  found  to  be  msuiliciently  accurate.  If  tire  hierarchy  of  models  is  chosen  in  such  a  way 
that  the  areas  where  the  lower  fidelity  models  fail  are  small  compared  with  the  size  of  the  design  space,  then 
the  procedure  described  above  should  be  quite  effective  in  producing  results  that  are  of  nearly  high-fidelity 
over  the  rnttir  design  space.  Our  experience  shows  that  this  is  the  ease  for  aerodynamic  performance:  the 
PASS  aerodynamic  module  is  quite  good  at  predicting  the  absolutely  host  wing  (lower  hound  estimate  mi 
the  Cq)  that  could  he  produced  if  considerable  design  work  were  done  on  ihe  configuration  (potentially 
using  adjoint  methods  and  a  liigh-diniensional  shape  parameterization),  However,  it  is  tumble  to  predict 
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Figure  !  1:  Database  of  PASS  (green),  A5G2/Pmiair  (red)  and  FE  (blue)  results. 


some  of  the  liner  details  of  aerodynamic  performance  anti  certainly  fails  when  transonic  effects  arc  present. 
Ao02/Panuir  is  also  unable  to  deal  with  transonic  flow  effects  but  produces  more  realistic  results  than  the 
PASS  analysis  as  the  actual  geometry  of  the  configuration  is  truly  accounted  for.  Finally,  the  Euler  models 
are  quite  good  predictors  of  the  aerodynamic  performance  of  the  complete  aircraft  as  long  as  viscous  effects 
are  not  dominant*  it  must  be  mentioned  that,  since  sonic  boom  1ms  not  been  considered  in  these  designs, 
the  Coarse  Euler  (CE)  evaluations  which  we  have  list'd  in  previous  work  would  be  sufficient  us  the  differences 
in  Cjj  between  CE  and  FE  calculations  were  found  to  be  small  (less  than  5  counts)  over  the  large  range  of 
variations  pursued  in  this  work. 


Second  Level  Optimization:  Complex  Configuration  Adjoint  Method 

Non-gradient  based  methods  (direct  searches,  simplex  method,  genetic:  algorithms)  have  been  shown  to 
work  with  a  wide  range  of  problem  types,  and  the  additional  complexity  required  to  compute  gradients, 
approximate  Hessians,  perform  line  searches  and  determine  the  optimal  steps,  is  not  needed.  These  properties 
have  a  tendency  to  make  this  search  procedures  simple  and  robust. 

In  addition,  these  search  procedures  are  able  to  handle  rather  noisy  and  large  design  spares,  milking 
them  a  reasonably  good  match  for  conceptual  design  procedures,  as  was  mentioned  in  the  Introduction. 
However  all  of  these  algorithms  tend  to  require  a  very  large  number  of  function  evaluations  for  convergence 
and  therefore  their  computational  cost  can  be  very  high. 

On  the  other  hand,  if  a  design  problem  shows  a  smooth  response  to  the  variation  in  the  design  variables 
and  gradient  information  is  readily  available  and  can  be  obtained  inexpensively  (as  is  the  case  with  the 
adjoint  method)  gradient-based  optimization  techniques  can  be  shown  to  Stave  significant  advantages  over 
non-gradient  search  procedures. 

In  our  approach  we  are  seeking  to  combine  the  advantages  of  both  gradient-  and  Tum-gmdiejit.- based  op¬ 
timization  procedures.  As  mentioned  earlier;  PASS  uses  a  simplex  method  and  is  able  to  produce  reasonable 
designs  (using  a  maximum  of  around  20-25  design  variables)  even  with  very  large  variations  of  the  design 
variables.  Once  the  simplex  method  has  converged  to  an  optimum  (local  or  global)  we  may  limit  ourselves  io 
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Cruise  Mach 

1.6 

Rang* 

um  inui 

BFL 

6.500  ft 

Minimum  static  margin 

0,0 

Alpha  limit 

15° 

MTOW 

$6,875  lbs 

Table  l:  Performance  requirements  for  optimized  baseline  configuration, 

smaller  changes  in  the  configuration.  These  changes  are  more  likely  to  result  on  well-behaved  design  span's 
that  can  be  tackled  with  an  adjoint  procedure  and  a  gradient- based,  optimization  algorithm. 

In  this  second  level  of  our  optimizations  we  limit  ourselves  to  modifications  in  the  twist  and  camber  of 
Hie  wing,  while  maintaining  the  same  wing  planform.  fuselage  and  relative  positioning  of  the  nacelles  and 
empennage. 

L  wo  different  tools  are  available  for  this  portion  of  the  overall  optimization; 


1,  SYN87-SB.  A  single-block,  wing-body  Killer  adjoint  optimization  code  that  uses  the  MPSOL  SQP 
algorithm  for  the  optimization  with  t>r  without  constraints.  SYN87-SB  allows  For  arbitrary  changes 
to  the  shape  of  the  fuselage  and  wing  and  is  able  to  enforce  thickness,  curvature,  and  fuel  volume 
constraints. 

2,  S YN 1 07- MB.  A  multi- block,  complete  configuration,  HANS  adjoint  optimization  code  that  also  uses 
the  NPSOL  SQP  algorithm  for  optimization  and  that  allows  similar  geometry  controls,  cost  functions, 
and  constraints  as  SYN87-SB,  but  that  can  be  made  to  treat  arbitrarily  complex  geometries  such  as 
the  complete  aircraft  configurations  that  are  the  subject  of  this  work. 

Optimization  Results 


For  subsequent  design  work,  an  optimized  baseline  geometry  was  generated  by  running  the  standard  version 
of  PASS  for  a  mission  with  the  performance  objectives  summarized  in  Table  1  Mission  requirement}!  mid 
geometric  constraints  for  the  baseline  configuration  were  based  on  numbers  that  were  felt  to  be  representative 
of  current  industry  interest.  The  value  of  the  MTOVV  is  the  result  of  the  optimization  as  this  was  the  objective 
function  of  the  design.  As  mentioned  before,  in  an  effort  to  generate  an  aircraft  achievable  using  current 
levels  of  technology,  advanced  technology  assumptions  were  kept  to  a  minimum. 

The  values  of  the  design  variables  for  the  resulting  baseline  configuration  (which  are  also  used  as  starting 
points  for  subsequent  designs)  arc  provided  in  Table  2,  Note  that  the  values  highlighted  in  red  were  not 
allowed  to  vary  during  this  initial  optimization.  In  addition  to  these  variables,  G  variables  representing  the 
radii  of  fuselage  stations  located  at  5%,  10%,  15%.  62.5%.  75%,  and  87.5%  of  the  fuselage  length  were  added 
to  allow  for  performance  improvements  aiid  to  maintain  cabin  and  cockpit  compartment  constraints.  Finally 
wing  section  changes  were  allowed  at  three  defining  stations.  The  twist  at  the  root/symmetry  plane  section, 
the  leading  edge  crank  section  and  the  tip  section  were  allowed  to  vary.  Furthermore,  the  value  of  the 
maximum  camber  and  the  location  of  maximum  camber  were  also  allowed  to  change  at  the  first  two  wing 
defining  stations.  This  makes  up  for  an  addition  7  design  variables  for  the  wing. 


High-Fidelity  Validation  of  Optimization  Results 

The  aerodynamic  performance  of  the  configurations  predicted  by  PASS  combined  with  the  response  surface 
fits  should  be  validated  with  our  high-fidelity  tool,  AirplanePlus.  As  mentioned  above,  although  the  results 
of  the  Filler  validation  are  slight  ly  different  from  what  the  lit  predicted,  the  optimized  configuration  shows  a 
good  improvement  in  aerodynamic  performance  while  satisfying  all  the  mission  requirements.  A  comparison 
of  surface  pressure  distributions  (for  both  the  lower  and  upper  surfaces  and  in  side  view)  is  shown  below  in 
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Wing  and  Tail  Geometry 


Wing  reference  area  (Sr<  /) 

l.  o7ft  rC 

Wing  aspect  ratio  (AR) 

4,0 

Wing  quarter-chord  sweep  (A) 

53.35° 

W  ing  i  ajier 

0,15 

Wing  dihedral 

3C 

1. ending  edge  extrnsi(  hi 

0.27ft 

1  Milling  edge  extension 

0.107 

Break  It  H  at  it  »n 

0,1 

Location  of  wing  root  LE 

0.2!)  i 

Root  section  t/e 

2.5% 

Break  section  t/r 

3.0% 

Tip  section  t/e 

2.5% 

Yri  i  ierd  i  ail  area  i  ,  S, ,  f  t 

0.125 

Vertical  (ail  Alt 

0.05 

Yerliral  (ail  A 

sc0 

Vert  n  ;d  rail  A 

0.0 

tb irixont a  1  tail  area  (  <  S  .  r  « 

0.0 

1  lori/i mt al  (ail  A If 

2.0 

ilurizrmlal  uni  \ 

5C° 

Mi iji/i anal  uni  \ 

0.3 

Fuselage  Geometry 


Maximum  fuselage  length 

125  ft. 

Minimum  cockpit  diameter 

GO  inches 

Minimum  cabin  diameter 

78  inches 

(  \d a n  leugi  ji 

25  ft 

Table  2:  Geometric  design  variables  for  design  optimization  anti  values  for  baseline  design. 
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Figure  15:  Summary  of  baseli  nr  con  figuration. 
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Baseline  configuration 


Opt imized  configural ion 


Figure  10:  Pressure  dist.rilmHon  plots  -  lower  surface, 


Figures  l(i,  17,  and  18  The  reader  should  notice  a  lightly  larger  fuselage  radius  around  the  nose  area  ami 
increased  wing  inhoard  sweep  which  has  reduced  the  shock  strength  significantly.  We  have  also  created  a 
multiblock  mesh  with  nearly  (1  million  nodes  for  the  optimized  configuration  that  we  intend  to  use  in  future 
work  for  complete  configuration  adjoint  designs.  Given  that  this  mesh  was  at  hand,  it  provided  us  with 
a  unique  opportunity  to  cross- validate  all  of  the  Euler  results  that  had  been,  up  to  then,  produced  with 
our  unstructured  flow  solver.  Airplane  Plus,  The  results  cannot  be  more  satisfying  since  they  provide  nearly 
identical  solutions  throughout  the  range  of  CLs.  This  is  important  because  the  airfoils  across  the  span  of  the 
configuration  have  rounded  leading  edges,  but,  because  of  the  low  thick noss-to-diord  ratios,  it  is  quite  hard 
to  put  enough  grid  resolution  around  the  leading  edge  using  nearly  isotropic  unstructured  meshes.  With  tile 
multiblock  approach,  anisotropic  cells  are  easily  created  around  the  leading  edge  and  can  resolve  the  effects 
oF  leading  edge  curvature  rather  nicely.  This  means  that  the  unstructured  Euler  solutions  arc  just  as  capable 
of  doing  so.  As  an  aside,  we  had  thought  earlier  on  that  some  of  the  discrepancies  between  the  Euler  solvers 
and  A502  were  due  to  the  inability  of  the  Euler  solver  to  capture  (with  a  coarse  leading  edge  mesh)  the 
leading  edge  suction.  This  drag  polar  comparison  seems  to  indicate  that  this  is  not  the  ease, 

A  detailed  description  of  this  research  is  given  in  the  thesis  of  Seougirn  Choi,  entitled  “Multi-Fidelity 
Multi-Disciplinary  Design  Optimization  of  Supersonic  Business  Jets", 
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Baseline  eogtfiguratbii  Optimized  configuration 

Figure  17:  Pressure  distribution  plots  -  upper  surface* 


O  pt  it 1 1  ized  cot i  I  i  guration 

Figure  18:  Pressure  (list  ritmt  ion  plots  -  side  view. 


22 


Algorithms  for  Automatic  Feedback 
Control  of  Aerodynamic  Flows 


An  airplane,  by  its  very  nature  of  design,  is  meant  to  bo  a  flow  control  device.  This  becomes  clear  when 
one  considers  a  steady  air  flow  with  and  wit  limit  the  airplane.  The  very  presence  of  the  airplane  alters  the 
How  pattern,  significantly  so.  Through  the  ages,  the  primary  goal  of  the  aerodyna midst  has  been  to  design 
airplanes  such  that  they  meet  certain  performance  criteria.  This  could  be,  for  example,  the  maximum  range 
of  the  airplane  or  the  drag  at  cruise  conditions.  Lower  drag  immediately  translates  to  lower  fuel  consumption 
and  hence  lower  operating  costs. 

The  aerodynamic  performance  of  an  airplane  is  determined  by  the  nature  of  the  surrounding  flow  Held 
under  given  flight  conditions.  The  most  important  motivation  for  flow  control  arises  thus: 

If  it  becomes  passible  to  control  the  nature  of  the  surrounding  fluid  flow ,  then  it  ts  conceivable,  that 
both  the  operating  envelope  and  the  aerodynamic  perfonrumce  of  an  airplane  within  that  rnvehju  ran  hr 
s  igr i  i f im  n  l  ty  e n  h  a n  ced . 

Exploring  techniques  for  active  flow  control  has  been  the  central  theme  of  this  research.  All  previous 
attempts  at  flow  control  have  either  involved  designing  simplistic  controls  for  complex  problems  or  complex 
frt  dhark  based  controls  fur  simple  problems.  Problems  like  separation  control,  drag  reduction  and  control  of 
t  ho  vortex  shedding  frequency  in  the  How  past,  a  cylinder  have  all  been  controlled  using  open  loop  controllers. 

Closed  loop  control  has  been  demonstrated  only  on  simplistic  models  derived  from  simulation  or  exper¬ 
iment. 

An  Ideal  Flow  Control  Law  should  have  the  following  properties: 

L  Broadly  applicable:  we  are  looking  for  an  algorithmic  framework  for  generating  flow  control  laws  for 
a  variety  of  problems.  The  development  of  such  a  fr  amework  would  enable  easy  analysis  and  design  of 
control  laws  for  a  variety  oi  flow  control  problems. 

2.  Scientific:  the  control  laws  should  be  based  on  a  realistic  model  of  the  fluid  system. 

3.  Robust  :sihould  account  for  variability  in  measurement,  actuation,  etc.  This  would  mean  that  the  control 
u  should  he  feedback  based 

«  -  F(*),  (4) 

where  x  is  rhe  current  system  state. 

Our  goal,  therefore,  is  to  develop  feedback  based  control  laws  that  are  derived  from  a  realistic  representation 
of  the  flow.  We  try  to  make  sure  that  the  framework  is  as  generic  as  possible.  lending  easy  extension  In  a 
variety  of  situations.  We  then  discuss  specific  applications  of  the  control  law  thus  derived,  including  control 
of  Flutter. 


Flow  Control  using  Adjoint  Sensitivities 

Following  along  the  lines  of  the  argument  made  in  the  previous  section,  the  Flow  Control  problem  can  he 
posed  as  the  following  optimization  problem: 

Minimize:  /(w,u)f  (5) 

where  w  is  the  vector  of  state  variables  and  u  is  the  vector  of  control  variables,  /  represents  the  quantity 
being  controlled.  It  could  lie  a  measure  of  the  drag,  expressed  as  a  difference-,  between  the  current  surface 


pressure  distribution  mid  an  ideal  surface  pressure  distribution  that  has  minimum  drag.  /  could  also  he  a 
measure  of  the  deviation  from  an  equilibrium  condition.  For  example,  in  the  case  of  flutter  control.  I  \> 
chosen  to  be  a  measure  of  the  plunge  and  pitch  of  the  wing.  Minimizing  /  in  this  context  is  equivalent  to 
bringing  the  system  back  to  equilibrium. 

The  state  vector  w  consists  of  the  Elder  state  variables  at  each  Finite  Volume  in  the  domain.  Thus,  if 
there  are  a  million  cells  in  the  domain*  the  dimension  of  w  is  five  million  for  a  3-D  Mow.  The  control  vector 
u  in  our  case,  consists  of  the  surface  normal  mass  fluxes  at  every  cell  along  the  surface  of  the  wing.  The 
dimensionality  of  the  control  vector  u  is  much  smaller  than  that  of  the  state  vector  w, 

The  optimal  control  u  is  one  that  minimizes  t he  cost  function  (5).  The  first  derivative  of  the  cost 
function  with  respect  to  the  control  variables  u  is 

dl  dl  1  Ow  dl 
fin  Ow  On  f  Ou 

While  the  second  term  is  fairly  straightforward  to  evaluate,  the  first  is  not.  This  is  because  the  state  w  and 
the  control  u  are  related  by  the  Euler  equations  which  are  of  the  form 

R[ w.  u)  =  0.  (7) 


In  the  Adjoint  framework,  we  do  not  evaluate  the  partial  derivative  ^  explicitly.  We  first  note  that 
Equation  (0)  is  identically  zero.  Therefore,  it  is  permissible  to  multiply  it  by  a  Lagrange  multiplier  t  and 
add  it  to  Equation  (0).  This  gives 
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Rearranging  the  terms  in  Equation  (8)  we  get 
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\ow  ll*  is  an  arbitrary  multiplier. 
Adjoint  equation 


We  can  choose  it  to  make  the  coelficienl  of  ^7  zero.  This  is  called  the 
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Ilie  expression  for  the  gradient  of  the  cost  function  then  becomes. 
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Virtual  Aerodynamic  Shaping 

We  lest  the  Mow  control  algorithms  formulated  in  the  previous  sections  on  a  Virtual  Aerodynamic  SHuping 
problem.  Here,  we  lake  an  aerodynamic  configuration,  and  try  to  make  it  behave  like  another  aerodynamic 
couliguralion  under  the  same  free  stream  conditions  by  implementing  surface  flow  control, 

A  basic  review  of  ideal  fluid  aerodynamics  reveals  that  including  a  mass  source  on  the  surface  of  an 
air  foil  has  the  effect  ofitu  reusing  the  curvature  and  including  a  mass  sink  has  the  opposite  effect.  Given  that, 
for  a  pre-determined  performance  measure,  every  operating  condition  has  an  optimum  shape,  it  follows  tlmt 


a  combination  of  a  shape  that  is  optimal  for  one  operating  condition,  combined  with  suitably  placed  sources 
and  sinks  will  result  in  a  system  that  meets  optimum  performance  criteria  for  a  wide  range  of  operating 
conditions. 

For  the  purposes  of  this  study,  steady  jets  are  used  in  order  to  simplify  numerical  modeling  Mass  (lows 
are  prescribed  at  the  wall,  and  the  jets  are  modeled  so  as  to  satisfy  the  normal  velocity/flux  conditions  at 
the  wall.  In  addition,  the  nett.  mass  flow  through  the  wall  is  assumed  to  he  zero. 


(Id) 


Feedback  Nature  of  Adjoint  Based  Control 


The  Fe.vdhar.k  nature  of  the  Adjoint  based  control  laws  thus  derived  become  dear,  when  the  Adjoint  boundary 
condit ions  are  examined. 


The  Adjoint  Boundary  Conditions  for  Virtual  Aerodynamic  Shaping 


We  postulate  that  I  lie  behavior  of  any  aerodynamic  configuration,  at  any  freest  ream  condition,  is  defined 
by  its  surface  pressure  distribution  Thus,  when  we  talk  about  virtually  shaping  one  configuration 

to  behave  like  another,  what  we  realty  mean  to  do  is  to  places  mass  sources  and  sinks  at  suitable  locations 
on  1  he  configuration  such  that  the  resulting  pressure  distribution  P((,  p)  is  the  same  as  the  target  pressure 
dist  ribution. 


To  achieve  tins,  we  try  to  minimize  a  cost  function  of  the  form 

i~~  I  (p- Pi)2<nh ■ 

1  M 


(14) 


whine  F  is  the  pressure  at  the  surface  with  no  controls,  and  Pa  is  the  desired  target  pressure.  Observe  that 
this  cost  function  is  a  special  instance  of  the  more  generic  cost  function  described  in  Equation  (5).  We  can 
choose  our  computational  co-ordinates  such  that  tq  =  —  0,  and  the  normal  direction  is  along  cV  When 

the  blowing  velocities  prescribed  at  the  wall  are  very  small,  the  Adjoint  boundary  condition  reduces  to 
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For  Virtual  Aeivdynamiv  Shaping,  this  is  simply 


(  P  -  Pa  )  -  V'2^21  +  A- 1/'4 j$23  «  (16) 

t  he  corresponding  Adjoint  gradient  can  then  finally  he  expressed  implicitly  by  the  equation 

O  =  -  (\/ji  +  '^.'w  +  V’:j  v  +  V'4'W  +  V-'fj  (  E  +  )  )  <''(>‘h\<!8t  ■  (IT) 

It  can  lie  seen  from  Equation  ( 1C)  that  the  Adjoint  boundary  condition  is  dependent  on  the  difference  between 
the  desired  state  of  the  system  J\i  and  t.he  c  urrent  state  of  the  system  /\  Moreover,  Equation  (17)  clearly 
shows  that  the  Adjoint  (dadienl  depends  only  on  the  flow  variables  at  the  boundaries.  Thus,  it  is  clear  that 
the  Adjoint  based  control  thus  derived  is  Feedback  based. 

It  should  he  noted  however,  that  the  Adjoint  equation  is  solved  using  the  computed  state  of  the'  flow 
field  and  not  the  measured  state  Tins  bypasses  the  need  for  global  measurements  of  t  he  flow-field. 
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Results 


2-D  results 

Aii  RAK-82  airfoil  was  optimized  for  minimum  drag  at  a  Mach  number  of  0.77.  The  airfoil  was  constrained  to 
operate  at  a  Cl  of  0.6,  The  pressure  distribution  of  the  optimized  section  was  used  as  the  target  distribution 
for  the  flow  control  rase,  where  blowing  and  suction  is  used  to  mimic  the  shape  changes  that  lead  to  the 
desired  pressure  changes.  The  flow  calculations  were  done  on  a  102  x  22  grid. 

The  original  (solid)  and  optimized  (dotted)  airfoil  are  shown  in  Figure  10,  The  blowing  and  suction 


Figure  10:  RAE-82  optimized  for  minimum  drag  at  Mach  0.77:  Original  (solid)  and  Optimized  (dotted) 


velocities  that  produce  the  same  pressure  distribution  are  shown  in  Figure  20.  As  expected,  the  controller 
implements  blowing  on  the  lower  surface  and  suction  on  the  upper  surface.  This  corresponds  very  well  with 
the  shape  change  being  represented,  where  the  curvature  is  increased  on  the  lower  surface  and  decreased  on 
l  he  upper  surface. 


Figure  20:  RAE-82:  Flow  control  velocities  on  the  lower  (left)  and  upper  (right)  surfaces  for  Virtual  Arm 
dyntntrir  Shaping 


Hie  Pressure  distributions  before  and  after  applying  flow  control  are  shown  in  Figure  21  It  can  be  seen 
that  l  he  How  control  algorithm  derived  for  the  virtual  aerodynamic  design  ease  achieves  the  desired  results. 


Figure  21:  Pressure  distributions:  target  (solid)  and  actual  (dotted)  before  (left)  and  after  {rigla  )  How 
control 


Reduction  in  the  number  of  actuators:  Design  Trade-off 

The  results  included  in  the  previous  section  were  for  the  case  where  blowing  mid  suction  is  implemented 
continuously  along  the  surface  of  the  airfoil.  Implementing  this  is  not  practical.  We  therefore  try  to  reduce 
the  number  of  actuators. 

We  would  prefer  to  automatically  arrive  at  t  he  optimum  number  of  actuators  needed,  and  their  locat  ions. 
We  do  this  by  looking  at.  the  Adjoint  gradient.  The  Adjoint  gradient  represents  the  sensitivity  of  the  cost 
function  with  respect  to  the  control  variables.  The  numerical  values  of  the  gradient  derived  thus  indicate 
which  controller  locations  are  most  effective  and  which  controller  locations  are  least  effective. 

In  order  to  frame  this  mathematically,  we  chose  to  include  all  locations  where  the  control  input  required 
was  ai  least  70  percent  of  that  when'  l-he  effort  was  maximum,  and  set  the  blowing/ suet  ion  velocities  at  all 
other  locations  to  zero. 

if  (t)  <  0-7  max  /«/„  ,  /*/„ (0=0.  (18) 

It  can  be  seen  from  Figures  22  and  22  that,  suction  control  is  applied  only  between  about  5  percent  chord 
and  30  percent  chord  on  the  upper  surface,  and  no  control  is  applied  otherwise.  The  magnitude  of  suction 
required  is  about  the  same  as  that  hi  the  continuous  control  case  (20). 


Design  TV  ado  off:  The  results  of  this  experiment  can  be  seen  in  Figure  23.  It  can  be  seen  that  the  desired 
pressure  distribution  is  almost  obtained.  The  match  between  the  desired  and  actual  pressure  distributions 
are  quite  close  especially  at  the  leading  edge  where  the  control  is  applied.  In  general,  we  make  a  compromise 
when  we  move  from  infinite  dimensional  control  to  limit*  dimensional  control.  The  l radc-ofT  between  i  br¬ 
ieve]  of  control  accuracy  desired  and  the  munber/location/size  of  the  controllers  that  can  he  implemented 
is  a  design  choice.  The  current  section  merely  presents  the  algorithm  that  one  would  use  to  study  the 
ramifications  of  such  a  tradeoff. 

3-D  Results 

f  inally,  we  check  to  see  if  we  can  achieve  similar  results  in  3  dimensions.  Pile  Surface  Pressure  distribution 
of  an  ONKRA  M(>  wing,  constrained  to  operate  at  a  O.  of  0,3  and  a  Mach  number  of  0.81  is  shown  in  Figure 
2L  We  try  to  achieve  the  same  surface  pressure  distribution  aT  the  same  freest  ream  conditions  for  another 
wing  that  has  a  NACA  0012  section.  The  flow  calculations  are  performed  on  a  192  x  32  x  18  grid.  3t  can 
be  seen  from  Figure  21  that,  after  5  control  iterations,  the  surface  pressure  distribution  resembles  that  of 
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1  1 


Figure  22:  RAE-82:  Flow  control  velocities  oil  the  lower  (left)  and  upper  (right)  surface  for  Virtual  Atroth/- 
namic  Shaping  reduced  number  of  actuators 


Figure  23:  Pressure  distrihut ions:  forget  (solid)  and  actual  (dotted)  after  How  control  -  Reduced  number  of 
Actuators 
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the  ON  ERA  MG.  The  pressure  distributions  along  the  front  portions  of  tho  wing  an*  almost  identical  The 
pressure  distributions  along  the  trailing  edge  exhibit  a  slight  difference.  The  original  pressure  distribution 
on  the  NACA  0012  wing  is  shown  in  the  dotted  linos. 
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Figure  24:  Cp  distribution  over  the  surface  of  an  Onera  MG  wing  at  M  ~  0,84  and  —  0  3  (left)  and 
Virtual  Aerodynamic.  Shaping  of  a  NACA  0012  wing  to  match  the  surface  pressure  distribution  of  an  Onera 
MG  wing  at  M  =  0.84  and  (\  —  0.3  (right) 


An  Algorithmic  Approach  to  Flutter  Control 

The  structural  design  of  an  airplane  is  guided  by  static  and  dynamic  factors.  The  more  stringent  constraints 
on  the  structural  design  are  due  to  dynamic  loads,  caused  by  aero-elastic  interactions.  One  of  the  most 
commonly  encountered  problems  in  aeroelasticity  is  flutter  [  1  ] ,  a  term  that  is  used  to  recognize  the  transfer 
of  energy  from  unsteady  aerodynamics  associated  with  the  surrounding  fluid  to  the  wing  structure,  resulting 
in  rapidly  divergent  behavior.  If  flutter  can  be  controlled  at  cruise  speeds*  we  can  design  lighter  wings 
and  consequently  more  efficient  airplanes.  It  is  therefore,  in  the  aircraft  designer’s  best  interest  to  design 
innovative  ways  in  which  flutter  can  be  controlled  without  making  the  resulting  structure  too  heavy. 


0.1  8 

2-D  Flutter  Control  In  the  present  section  we  will  investigate  the  aeroelastir  behavior  and  control  of  a  2-D 
airfoil  whose  schematics  is  shown  in  Figure  25.  A  2-D  airfoil  model  can  be  shown  to  be  a  fair  represent  at  a  in 
for  flutter  prediction  as  shown  by  Theodorson  and  Garrik  [3]  of  a  straight  wing  of  a  large  span  by  giving  it 
the  geometric  and  inertial  properties  of  the  cross-section  three  quarters  of  the  way  from  the  centerline  to  the 
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winy;  tip.  The  equations  of  motion  of  this  simple  system  can  he  shown  to  be  as  follows. 


Figure  25:  Typical  Section  Wing  Model  Geometry 


rnh  +  5An  +  Kfi  h  —  —  L  (19) 

SJi  +  Itx&  +  Kao  -  (20) 

fWi  anti  Ka  are  representative  of  the  bending  and  torsional  si  illness  of  the  wing  about  its  elastic  axis. 

the  elastic  axis  is  the  Inens  of  points  about  which,  if  a  force  is  applied,  doesn't  result  in  any  rotation 
about  that  point,  rn  and  Itr  are  the  mass  and  moment  of  inertia  of  the  wing  section  about  the  elastic  axis. 

Sn  i>  the  coupling  term  which  depends  on  the  relative  position  of  the  center  of  gravity  and  the  elastic  axis 

We  assume  I  hat  the  structural  properties  arc  fixed  and  we  have  some  amount  of  control  of  the  right 
luind  sides  of  Equations  (19)  and  (20)  via  blowing  and  suction  The  objective  is  to  find  a  suitable  control 
law  which  will  modify  the  aerodynamic  terms  so  as  to  prevent  flutter. 

Computational  simulation 

!  hr  flow  is  simulated  by  solving  the  unsteady  Euler  equations.  The  Eider'  eq nations  are  solved  using  a  dual 
time  stepping  method,  using  a  third  order  backward  difference  formula  in  time,  and  a  symmetric  Gauss  Seidel 
scheme  for  solving  the  inner  iterations.  The  above  mentioned  flow  simulation  code  is  integrated  with  a  two 
degree  of  freedom  structural  model  given  by  Equations  (19)  and  (20).  The  coupled  aero-structural  system  is 
integrated  using  the  Nt wm m k  scheme.  The  simulation  techniques  are  discussed  in  detail  in  Chapter  2  and 
Appendix  B  of  the  thesis  of  Palau iappanV 


System  Linearization  and  Model  Order  Red  net  ion 

In  Equations  (19)  anti  (20),  the  structural  parameters  are  constant.  The  lift  L  and  the  moment  M  arc 
complex  nonlinear  functions  of  the  system  state  w.  o.  d.  h  and  h.  Moreover,  o,  ri,  h  and  h  are  itself 
functions  of  the  system  state  w.  Here  the  state  w  is  the  vector  consisting  of  all  the  Eufcr  stales  at  all  finite 
volumes  used  in  the  simulation.  Thus 
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“  />(w.h). 

(21) 

M 

-  A/(w,u). 

(22) 

Linearizing  about  the  nominal  operating  point. 

we  get 
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It  should  be  noted  i  hni  for  a  simulat  ion  wit h  one  million  finite  volumes,  t  lie  dimension  of  w  is  four  million  for 
a  2*1)  simulation  and  five  million  For  a  5- 1 )  simulation,  l  lnis  evaluating  the  above  derivatives  is  a  formidable 


computational  challenge.  ft  is  also  important  to  recognize  that  not  all  the  derivatives  me  significant  in  the 
above  representation.  Consider,  for  example,  a  cell  in  the  far-held.  The  value  o(  the  state  variable  tliere 
is  not  going  t o  change  by  much,  however  rapid  the  oscillations.  Therefore,  it  is  of  very  little  use  evaluating 
these  derivatives  in  our  linearized  model. 

Instead,  we  choose  to  obtain  a  suitable  reduced  order  model  that  captures  the  essential  physics.  The 
most  obvious  reduction  t  hat  we  can  obtain  is  in  terms  of  a,  a,  h  and  it,  We  therefore  work  with  a  mode!  of 
the  form: 


L 


M 


f  f  .  ,  ,  ,  ;  t)L‘ 

4-  /jqfl  4-  LiJi  4-  Li  h  +  — —  u, 

a  u 

(UfT 

4-  AL\t\  4-  ,\lf,h  4-  A l ill  +  — ■ —  u 

du 


(25) 
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Equations  (25)  and  (20)  assume  that  the  nominal  values  of  n,  d,  h  nucl.  h  and  u  are  zero,  respectively.  Tims 
for  l lie  flutter  control  problem  being  studied,  the  following  state  vector  is  used: 

x  =  [or  a  ft  h]7  (27) 


System  Identification:  Evaluation  of  Sensitivities 

In  our  aero-str  neutral  model  (19)  and  (20).  the  lift  L  and  the  moment  M  depend  on  the  complete  system 
state  w.  However,  using  a  full  order  state  model  to  design  a  controller  is  not  feasible,  given  the  extremely 
high  dimensionality  of  the  system.  We  therefore,  formulate  a  reduced  order  model  of  the  system  as  shown  in 
Equations  (25)  and  (2(i).  Tn  order  for  tins  model  to  be  complete,  we  need  to  evaluate  the  sensitivities  with 
respect  to  the  reduced  order  state  x  and  the  control  variables  u. 


Sensitivities  with  respect  to  the  state  variables 

The  sensitivities  of  the  lift  and  moment  with  respect  to  the  state  variables  are  evaluated  in  two  different 
ways. 


Tbeodorsen  theory:  First,  we  use  theoretical  results  horn  Theodorsen  [I  ] .  Theodorsen  theory  assumes 
that  the  airfoil  under  consideration  is  thin,  and  is  oscillating  in  an  incompressible  flow.  Under  those  consid¬ 
erations 


L„  =  npv^c  . 

,  irpv oo/r2 

h,  =  (J  . 

,  Lfl  =  ■n(>vQO(\ 

■  M&  =  U, 

Mh  =  0  , 

, ,  irjM'ocf2 

■  Mk  =  . 

Hero  {>  is  the  freest  ream  density,  v ^  is  the  freest  ream  velocity  and  c  is  the  chord  of  the  airfoil. 


Least-Squares  Method:  In  the  second  method,  we  evaluate  the  sensitivities,  by  studying  the  unforced 
response  of  a  pitching  airfoil,  and  then  estimating  the  sensitivities  by  a  least-squares  technique.  The  aero- 
structural  response  of  the  system  over  a  period  of  time  is  similar  to  the  unforced  response  reproduced  in 


Figures  27,  27,  28  and  28.  These  simulation#  provide  numerical  values  For 


n  =  /jf/) 

d  =  m 

h  =  f2(l) 

u  =  m 

l  =  Mt) 

M  -  Mt] 

We  now  try  to  fit.  the  data  thus  obtained  to  functions  of  the  form 


L  —  4  4-  Istill  4  I'jJ*  y 

A/  —  Aifi  ft  d*  A  /ri  (V  4  A  //i  /i  +  A  /  ^  h , 

Our  goal  is  to  evaluate  the  sensitivities  LCI ,  Att.  A/„.  A  A,.  .U/,  and  A/,(.  We  do  this  using  a  least- 

squares  technique. 

It  can  he  seen  from  the  simulation  results  that  both  techniques  work  quite  well  The  system  identification 
by  the  least -squares  technique,  works  slightly  better,  in  the  sense,  it  achieves  faster  stabilization.  This  r  an 
be  attributed  to  the  fact  that  this  represents  the  nonlinear  system  it  tore  closely. 


Sensitivities  with  respect  to  the  control  variables 

We  also  need  to  evaluate  the  sensitivities  of  L  and  A/  with  respect  to  the  blowing  and  suction  velocities  u. 
£  auc!  4  respectively. 

We  cio  tins  arc  using  an  Adjoint  method  as  outlined  in  Chapter  . 


Flutter  Control:  Formulation  of  the  Objective  Function 

We  can  define  1  he  flatter  velocity  as  that  point  where  wo  have  sustained  oscillations  of  the  system.  Let  as 
define  the  state  vector  x  as  follows 

X  =  [a  &  h  li]T  (28) 

The  control  vector  u  is  the  vector  of  blowing/suction  velocities  at  the  wall  The  dynamics  of  the  system  is 
represented  by  (19.)  and  (20).  For  the  purposes  of  designing  a  controller,  we  model  the  lift  L  and  the  moment 
M  using  a  reduced  order  model  as  presented  in  Equations  (25)  and  (26).  The  system  model  used  to  design 
a  controller  is  then 

mh  +  Sna  4  h'tji 

Snh  4  4  /C>o 

This  can  be  re-phrased  in  state  space  form  as  follows: 

A/x  =  /ix  4  Bn.  (29) 

Here  the  matrix  //  represents  the  sensitivities  of  the  state  vectors  with  respect  to  the  control  variables.  Tins 
can  be  obtained  by  solving  the  Adjoint  equations.  Inverting  A/,  we  get  a  system  of  the  form 

x  —  /lx  4  Bn.  (;j<)) 
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It  is  possible  to  design  a  controller  for  the  system  (30)  using  LQR  ta'i-liiiiqiMW  [2),  The  objective  of  the 
problem  is  to  control  tile  system  given  by  (.'ill),  so  that  the  final  value  of  the  state  vector  is  given  by 

xr  =  [of  l)  hr  0]T  (111) 


If  this  is  rephrased  as  an  optimization  problem,  the  objective  would  be  to  minimize  the  following  funrti 


Xf )?  Q(x  Xf )  4-  II J  Ru)  dt 


m 


where  Q  is  n  positive  semi-definite  weighting  matrix  and  R  is  a  positive  definite  matrix.  In  our  ease, 


Q  -  /. 

R  =  sli 


where  /  is  the  identity  matrix,  and  f  is  a  small  positive  constant.  R  is  required  to  be  positive  definite,  to 
ensure  that  the  control  computed  is  not  of  unreasonable  magnitudes. 


Backsubstitution  of  the  Control  Law  into  the  Nonlinear  System 

A  Feedback  control  gain  matrix  of  the  Stanford  LQR  form  is  then  derived  for  l  lie  flutter  control  problem.  Now, 
the  aero- struct u rnl  system  is  simulated  with  blowing  and  suction  control  applied  at  the  actuator  locations. 
1  he  magnitude  of  control  requited  at  each  actuator  location  is  given  by  the  control  gain  matrix  Kss 

u  =  (33) 

It  can  be  seen  that  this  control  law  was  successful  in  stabilizing  the  system.  The  results  are  presented  in  the 
next  section. 


Results 


The  following  experiments  were  conducted  on  a  symmetric  NACA  0012  airfoil  at.  a  freestream  Mach  number 
of  0.3.  A  1  fill  x  32  grid  was  used  for  the  CFD  simulation. 

The  structural  properties  were  chosen  as  follows:  ilk  -  GO,  M  —  GO,  Kft  ~  GO,  Kn  =  GO,  and  Sn  —  31), 
Our  nominal  rest  point  is  a  =0°  and  h  —  0, 

Adjoint  Gradients 

As  discussed  in  the  last  section,  the  Adjoint  method  is  used  to  find  Hie  gradients  of  lift  and  moment  with 
respect  to  the  control  variables,  namely  the  blowing  and  suction  velocities  cm  the  surface.  It  should  be  noted 
that  lids  is  done  using  a  steady  flow  assumption  about  the  nominal  rest  point  of  the  system.  We  used  a 
symmetric  NACA  0012  sect  Urn.  So  for  our  case,  this  nominal  rest  point  was  at  a  =  0.  and  h  —  0.  These 
gradients  are  shown  in  Figure  20, 


Application  of  Feedback  Control  to  the  Nonlinear  Flutter  Problem 

The  uncontrolled  and  controlled  aero-structural  simulations  are  represented  in  Figures  27,  and  28.  It  should 
be  noted  that  even  though  the  feedback  law  is  derived  from  a  linearized  model  of  the  system,  the  control  is 
applied  to  a  complete  nonlineai  model.  Two  different  methods  are  used  to  find  the  aerodynamic  derivatives. 
It  can  he  seen  that  the  least-squares  method  does  a  better  job  than  the  Theodorsen  method  for  Hut  lei 
control.  This  is  obvious  because  this  represents  the  nonlinear  system  more  accurately.  The  corresponding 
blow  ing/s  tie  Hon  velocities  are  shown  in  Figure  29,  It  should  be  noted  that  the  freestream  value  of  pqn  in 
our  simulation  was  1.  So  the  values  of  blowing  and  suction  controls  required  is  quite  smalt.  Moreover,  we 
need  zero  control  input  at  the  equilibrium  point,  which  is  what  we  desire. 


Figure  26:  Gradient  of  lift  (left)  and  moment  (right}  with  respect  la  control  mass  fluxes 
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Figure  27:  Variation  of  angle  of  attack  (degrees)  with  time:  controlled  and  uncontrolled  cases  (left)  ami 
Variation  of  Cm  with  time:  controlled  and  uncontrolled  cases  (right) 


Figure  2  H:  Variation  of  plunge  h/c  with  time:  controlled  and  uncontrolled  cases  (left)  and  Variation  off) 
with  time:  controlled  and  uncontrolled  cases  (right) 


Figure  29:  Blowing/ Sue  I  ion  mass  fluxes  at  the  Leading  Edge  (left)  and  "rime  step  refinement  studies  for  the 
variation  of  angle  of  attack  with  time  (right) 

Time  step  refinement  studies 

To  ensure  that  the  flutter  control  simulations  are  correct,  the  time  step  for  the  nonlinear  aero-structural 
solver  is  made  smaller  and  smaller  and  the  controlled  behaviour  is  observed.  It  can  be  seen  that  the  pattern 
of  variation  of  the  angle  of  attack  with  time  is  fairly  well  predicted  by  the  solver.  (See  Figure  29). 


Reduction  in  the  number  of  Actuators 


Our  next  step  is  to  specialize  the  control  law  thus  derived  to  work  when  the  number  of  ad  uafors  is  finite. 
It  was  found  that  flutter  could  be  controlled  with  as  few  as  four  actuators:  one  each  in  the  leading  and 
trailing  edges  and  one  each  in  the  middle  of  t lie  upper  and  lower  surfaces.  The  fact  that  there  are  only  four 
actuation  points  is  represented  by  zeroing  out  the  gradient  shown  in  Figures  20  and  2G  everywhere  except 
at  these  four  locations.  (Every  location  is  represented  by  a  small  cluster  of  CFO  cells  to  prevent  numerical 
instability  and  damping  of  the  actuation  values.) 

The  entire  procedure  outlined  in  the  previous  section  is  then  repeated  to  derive  the  feedback  gain  matrix 
/w  It  can  be  seen  from  Figures  30.  30,  31  and  31  that  the  matrix  has  non-zero  values  only  at  the  desired 
locations  of  the  controllers.  Consequently,  actuation  is  performed  only  at  these  sites.  This  is  equivalent  to 
eont  rolling  the  problem  with  a  finite  number  of  actuators. 

It  can  be  seen  from  Figure  32  that  flutter  is  controlled  successfully  even  wit  h  a  finite  number  of  actuators. 
This  is  an  important  result,  as  it  implies  that  this  system  can  la*  implemented  mi  a  practical  aerodynamic 
coufigurat  ion. 


3-D  Results 

We  t  hen  try  to  control  the  flutter  of  a  realistic  airplane  wing.  The  wing  is  i inswept,  and  the  cross-section  is 
that  of  a  0  percent  thick  airfoil  obtained  by  scaling  down  a  NACA  111)12  airfoil.  The  semi-span  of  the  wing 
is  11.5  inches,  and  t  he  chord  is  1.50  inches.  This  corresponds  to  an  aspect  ratio  of  about.  5. 

Structurally,  the  wing  is  modeled  as  a  plate  of  thickness  0,005  indies  that  is  placed  along  the  centerline 
of  the  wing-section.  The  density  of  the  material  of  the  wing  is  0.003*108  slug/sq.  inch.  The  Young’s  modulus 
is  9,818  x  IQ6  slug/ sq,  inch  and  the  torsional  rigidity  is  3.039  x  U)i}  slug/sq.  inch. 


Figure  3l>:  Coefficient  c>E  rotation  angle  vs,  actuator  number  in  the  feedback  gain  matrix  (left)  uml  Coefficient 
of  plunge  vs.  actuator  number  in  the  feedback  gain  matrix  (right) 


Figure  31:  Coefficient  of  rotation  angle  rate  vs,  actuator  number  in  the  feedback  gain  matrix  (left)  and 
Coefficient  ot  plunge  velocity  vs.  actuator  number  in  the  feedback  gain  matrix  (right) 


Figure  32:  Variation  of  angle  of  attack  (degrees)  with  time:  with  1  actuators 
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The  wing  was  operated  under  a  freest  ream  Mach  number  of  0.79  and  a  freestream  dynamic  pressure  of 
5211  Pascal. 

The  structure  is  modeled  using  50  plate  elements.  The  aerodynamic  simulation  is  dune  on  ri  96x32  *  In 
grid.  1(  can  be  seen  from  Figure  33  that  the  uncontrolled  system  diverges  fairly  rapidly.  In  the  time  frame 
considered,  the  plunge  diverges  from  a  negligible  amount  to  It)  percent  chord  very  quickly. 

Our  task,  now.  is  to  design  a  controller  using  the  techniques  described  in  the  previous  sections,  It  1ms 
been  shown  that  the  flutter  of  a  wing  can  he  studied  by  studying  the  dynamics  of  a  section  three  quarters 
of  the  distance  from  the  wing  center-line  to  the  tip.  We  identify  the  structural  properties  of  the  section 
located  at  this  point,  and  mode]  it  using  the  typical  section  wing  model,  discussed  previously.  Following  the 
techniques  in  the  previous  section,  we  quickly  derive  the  feedback  gain  matrix  t\^  for  this  section. 

We  make  the  assumption  that  this  matrix  is  valid  throughout  the  wing.  This  is  a  valid  assumption*  as 
the  control  is  implemented  in  a  feedback  fashion.  The  t ip  is  expected  to  go  through  the  maximum  deflection, 
and  therefore  will  be  subject  to  the  maximum  amount  of  control.  (Since  the  control  is  proportional  in 
nature).  The  root  does  not  move  at  all.  and  thus  there  is  no  control  applied  at  the  root.  The  results  of 
fins  simulation  are  shown  in  Figure  33.  It  ran  be  seen  that  the  control  law  thus  derived  is  successful  hi 
controlling  Mutter.  Hie  mass  Maxes  at  an  actuator  location  at  the  tip.  along  the  trailing  edge  arc  shown  in 
Figure  33.  Again,  it  can  be  seen  that  the  mass  fluxes  required  for  control,  when  compared  to  the  freestrouin 
mass  flux  of  pqn  —  1  are  very  small. 


Figure  33:  Variation  of  plunge  h/c  with  time:  controlled  and  uncontrolled  cases  (left)  and  Blowing/Snetimi 
mass  fluxes  at  a  trailing  edge  point  (right) 

In  the  controlled  case  the  wing  settles  into  n  limit  cycle  oscillation  of  small  magnitude  as  can  be  seen 
from  Figure  33.  This  is  in  spite  of  the  fact  that  an  approximate  structural  model  was  used  in  the  calculation 
of  the  control  law. 

This  research  is  described  in  detail  in  the  thesis  of  Kartlrik  Palanmppau,  entitle  "Algorithms  for  Auto¬ 
matic  Feedback  Control  of  Aerodynamic  Flows”. 
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